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Abstract 

The equations of general relativity in the form of timelike and null 
geodesies that describe motion of test particles and photons in Kerr space- 
time are solved exactly including the contribution from the cosmological 
constant. We then perform a systematic application of the exact solu- 
tions obtained to the following cases. The exact solutions derived for null, 
spherical, polar and non-polar orbits are applied for the calculation of 
frame dragging (Lense-Thirring effect) for the orbit of a photon around 
the galactic centre, assuming that the latter is a Kerr black hole for var- 
ious values of the Kerr parameter including those supported by recent 
observations. Unbound null polar orbits are investigated, and an analyt- 
ical expression for the deviation angle of a polar photon orbit from the 
gravitational Kerr field is derived. In addition, we present the exact solu- 
tion for timelike and null equatorial orbits. In the former case, we derive 
an analytical expression for the precession of the point of closest approach 
(perihelion, periastron) for the orbit of a test particle around a rotating 
mass whose surrounding curved spacetime geometry is described by the 
Kerr field. In the latter case, we calculate an exact expression for the 
deflection angle for a light ray in the gravitational field of a rotating mass 
(the Kerr field). We apply this calculation for the bending of light from 
the gravitational field of the galactic centre, for various values of the Kerr 
parameter, and the impact factor. 
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tMIFP-05-16, July 2005 



1 



1 Introduction 
1.1 Motivation 



Most of the celestial bodies deviate very little from spherical symmetry, and 
the Schwarzschild spacetime is an appropriate approximation for their gravi- 
tational field [2]. However, for some astrophysical bodies the rotation of the 
mass distribution cannot be neglected. A more general spacetime solution of 
the gravitational field equations should take this property into account. In this 
respect, the Kerr solution [3] represents, the curved spacetime geometry sur- 
rounding a rotating mass [4]. Moreover, the above solution is also important for 
probing the strong field regime of general relativity [5]. This is significant, since 
general relativity has triumphed in large-scale cosmology [8, 7, 9, 10], and in 
predicting solar system effects on planetary orbits like the perihelion precession 
of Mercury with a very high precision [1, 6] ^. 

As was discussed in [11], the investigation of spacetime structures near strong 
gravitational sources, like neutron stars or candidate black hole (BH) systems is 
of paramount importance for testing the predictions of the theory in the strong 
field regime. The study of geodesies are crucial in this respect, in providing 
information of the structure of spacetime in the strong field limit. 

The study of the geodesies from the Kerr metric are additionally motivated 
by recent observational evidence of stellar orbits around the galactic centre, 
which indicates that the spacetime surrounding the Sgr A* radio source, which is 
believed to be a supermassive black hole of 3.6 million solar masses, is described 
by the Kerr solution rather than the Schwarzschild solution, with the Kerr 
parameter [12] 

^^=0.52 (±0.1, ±0.08, ±0.08) (1) 

where the reported high-resolution infrared observations of Sgr A* revealed 'riui- 
escent' emission and several flares. This is half the maximum value for a Kerr 
black hole [13]. In the above equation J ^ denotes the angular momentum of 
the black hole (The error estimates here the uncertainties in the period, black 
hole mass (Mbh) and distance to the galactic centre, respectively; G is the 
gravitational constant and c the velocity of light.) 

Taking into account the cosmological constant A contribution, the general- 
ization of the Kerr solution is described by the Kerr -de Sitter metric element 
which in Boyer-Lindquist (BL) coordinates ^ is given by [18, 19]: 

^See also about the BepiColombo science mission on Mercury 
http: / /sci.esa.int/home/bepicolombo/ 

^J = ca where a is the Kerr parameter. The interpretation of ca as the angular momentum 

per unit mass was first given by Beyer and Price [15]. In fact, by comparing with the Lense- 

Thirring calculations [16] they determined the Kerr parameter to be: a = — , where Q 
and / denote the angular velocity and radius of the rotating sphere. 

^ These coordinates have the advantage that reduce to the Schwarzschild solution with a 
cosmological constant in the limit a — » 0, see [14]. 
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{cdt — asin^ 





\— (acdt - (r^ + a'^)d4>f 



(2) 



where 



A, 



A. 




(3) 



In a recent paper [11], we derived the timelike geodesic equations in Kerr 
spacetime with a cosmological constant by solving the Hamilton- Jacobi partial 
differential equation by separation of variables. Subsequently, we solved exactly 
the corresponding differential equations for an interesting class of possible types 
of motion for a test particle in Kerr and Kerr-(anti) de Sitter spacctimos. The 
exact solution of non-spherical geodesies was obtained by using the transforma- 
tion theory of elliptic functions. 

The exact solutions of the timelike geodesic equations obtained in [11] were 
applied to the following situations: 

Frame dragging from rotating gravitational mass. An essential property of 
the geodesies in Schwarzschild spacetime is that although the orbit precesses 
relativistically it remains in the same plane; the Kerr rotation adds longitudi- 
nal dragging to this precession. For instance, in the spherical polar orbits we 
discussed in [11], (where the particle traverses all latitudes, passes through the 
symmetry axis z, infinitely many times) the angle of longitude increases after a 
complete oscillation in latitude. This phenomenon, is in accordance with Mack's 
principle. 

More specifically, in [11] we calculated the dragging of inertial frames in the 
following situations, (a) Dragging of a satellite's spherical polar orbit in the 
gravitational field of Earth assuming Kerr geometry, using as radii, the semi- 
major axis of the polar orbit of the GP-B mission launched in April 2004. 
(b) Dragging of a stellar, spherical polar orbit, in the gravitational field of a 
rotating galactic black hole. 

It is the purpose of this paper to extend the analysis and applications of 
the exact solutions obtained in [11] to an interesting class of possible types 
of motion for a test particle in Kerr and Kcrr-(anti) de Sitter spacctimcs, as 
well as to derive the exact solutions of null geodesies in the same spacetimes 
and explore their physical implications. In the latter case, we apply the exact 
solutions obtained, to the following situations: 

*http://einstein.stanford. See also [17]. 
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(a) Dragging of a photon's spherical polar and non-polar orbit in the gravi- 
tational field of a rotating galactic centre black hole. 

(b) The deflection angle of a light ray from the gravitational fleld of a rotating 
black hole, for various values of the Kerr parameter and the impact factor. 

The material of this paper is organized as follows. In section 2 we review 
the derivation of the relevant geodesic equations. In section 2.1 we discuss 
the definition of the Lauricella's hypergeometric fmiction of many variables, as 
well as the integral representations that it admits, which are important in our 
exact treatment of geodesic equations that describe motion of a test particle 
and photon in Kerr-(anti) de Sitter spacetime. In sections 3.2,3.1 and 6, we 
solve exactly spherical polar or non-polar null geodesies with and without the 
cosmological constant. In the case of spherical polar photonic geodesies and 
for a vanishing cosmological constant the exact solution for the orbit is given 
by the Weierstrafi elliptic function. The exact expression for frame dragging 
is proportional to the real half-period of the Weierstrai3 modular Jacobi form. 
In the case of non-polar spherical photonic orbits, the exact expression for the 
Lense-Thirring precession of the photon is given in terms of a hypergeometric 
function of one variable and Appell's flrst generalized hypergeometric function 
of two variables Fi [20] . Assuming a vanishing cosmological constant and that 
the galactic centre is a supermassive rotating black hole, we apply the exact so- 
lutions obtained for the determination of Lense-Thirring precession for a photon 
in spherical polar and non-polar orbits around the galactic centre. The corre- 
sponding exact expressions in the presence of the cosmological constant are also 
derived and discussed in sections 3.2 and 7. 

In section 4 we perform a precise calculation for the deflection angle of a 
photon's non — spherical polar orbit from the gravitational Kerr field. In this 
novel case, the exact expressions obtained were written in terms of Lauricella's 
hypergeometric function Fjj. 

Timelike spherical polar (with A ^ 0) and non-polar (with A = 0) orbits are 
treated in sections 5, 8 respectively. 

In sections 8.4 and 8.5 we study the exact solution of non-spherical timelike 
and null equatorial orbits respectively. The amount of relativistic precession for 
a test particle in a timelike orbit, confined to the equatorial plane, in the presence 
of rotation of the central mass is given in terms of Appell's first hypergeometric 
function of two variables Fi . On the other hand, the exact expression for the 
deflection angle of a photonic orbit from the Kerr gravitational field surrounding 
a rotating central mass is given in terms of Appell's Fi hypergeometric function 
and Lauricella's fourth hypergeometric function Fjj of three variables [21]. We 
use section 9 for our conclusions. In the appendices, we collect some of our 
formal calculations, as well as some useful properties of Appell's hypergeometric 
function and definitions of genus-2 theta functions. 
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2 Separability of Hamilton-Jacobi's differential 
equation in Kerr-(anti) de Sitter metric and 
derivation of geodesies. 

In the presence of the cosmological constant it was proved in [11] the important 
result that the Hamilton- Jacobi differential equation can be solved by separation 
of variables. Thus in this case, the characteristic function separates and takes 
the form [11] 



/ 



W = -Ect + L(f)+ d6 



'S2 [(r2 + a^)E - aLY - Ar{i^\^ + Q + E^{L - aEf) 
+ / ^ dr 



By differentiating now with respect to constants of integration, Q,L,E,iJ, 
we obtain the following set of geodesic differential equations 



r dr _ r dO 

J J ~7W 



2 rft _ E\r^+a^) [jr^ + a^)E - aL] _ oE? {aE sm^ 9 - L) 
dX ~ Ar Ae 

where 

R' := E?[{r^ + a'')E-aL\' -Ar{ii''r^ + Q + E?{L-aEf) 

e' := [g + (L - aEfE? - /x^a^ cos^ 9] Ae - E? '^""'f ~ (5) 

sin 9 

The first line of Eq.(4) is a differential equation that relates a hyperelliptic 
Abelian integral to an elliptic integral which is the generalisation of the theory 
of transformation of elliptic functions discussed in [11], in the case of non-zero 
cosmological constant. The mathematical treatment of such a relationship was 
first discussed by Abel in [34]. 

Assuming a zero cosmological constant, as was shown by Carter, one gets 

W = -Ect + J ^dr + J Ved9 + Lcp (6) 
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where 



e :=Q 



a2(^2 - E^) + 



sm^e 



cos^O 



and 



with A 



R — + a^)E -aLf-A [fx^r^ + {L - aEf + Q] 



(7) 



(8) 



2GMr 



Also L are constants of integration associated with 
the isometries of the Kerr metric. Carter's constant of integration is denoted 
by Q. Differentiation of (6), with respect to the integration constants E, L, Q, fj. 
leads to the following set of first-order equations of motion [37] : 



where 



2 cdt 
2 dr 

fTx 



A 

±\/e 



—P - aE 
A 



-P-a{aE sin^ 6 - L) 



sin^6» 



P := E{r^ +0^) - ah 



(9) 
(10) 



Null- geodesies are derived by setting = 0. 



2.1 Lauricella's multivciriable hypergeometric functions 

Giuseppe Lauricella, building on the work of Appell who had developed hyper- 
geometric functions of two variables, investigated in a systematic way, multiple 
hypergeometric functions at the end of the nineteenth century [21]. He defined, 
four functions which are named after him and have both multiple series and 
integral representations. In particular, the fourth of these functions, denoted 
by Fy), admits integral representations of importance in our exact treatment of 
geodesic equations, which describe motion of a test particle in Kerr-(anti) de 
Sitter spacetime. 

The fourth Lauricella function of m-variables is given by 



oo 



n iTjriiH h"m Ujni ' ' ■ UJim 



E 



X;(a;,ni H n„)(/3i,ni) • • • (/3 



where 



/3 = (/?!,- (11) 
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The Pochhammer symbol {a)m = (a, m) is defined by 

^ r(a + m) _ f 1, if m = 0, 

r(a) ~ I a(a + l)---(a + m-l), if m = 1, 2, 3, • • • 

The series admits the following integral representation 

Fn{a,^,r,z) = ['e-\l-tr-"-\l-zrtr''^---{l-Zmt)-''-dt 

r(a)r(7 - a) Jo 

(12) 

which is vahd for Re(a) > 0, Re(7 — a) > 0. It converges absolutely inside the 
m-dimensional cuboid: 

\zj\<l, (j = l,---,m) (13) 

3 Spherical polar null geodesies 

Depending on whether or not the coordinate radius r is constant along a given 
geodesic, the corresponding particle orbit is characterized as spherical or non- 
sphcrical respectively. In this section, we will concentrate on spherical polar 
photon orbits with a vanishing cosmological constant. We should mention at 
this point the extreme black hole solutions a = 1 of spherical non-polar photon 
geodesies obtained in [22] in terms of formal integrals. 

The exact solution of the corresponding timelike orbits and their physical 
applications have been derived and investigated in [11]. 

Assuming a zero cosmological constant, r = rf, where rj is a constant value 
setting /X = and using the last two equations of (9) we obtain: 



del) _^-aE + L/ sin^ i 

where O now is given by 



(14) 



e = Q- [-a^E^ + cos^ e (15) 

sin 6 

It is convenient to introduce the parameters 

$ := L/E, Q := Q/E^ (16) 
Now by defining z := cos^ 9, the previous equation can be written as follows, 

1 faP $ 1 

2 ^z^a-z-^{a + (3) + Qz lA \-z] ^ ' 

where 

a:=-a2, /J := Q + $2 (^g) 
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It has been shown [23] that a necessary condition for an orbit to be polar 
(meaning to intersect the symmetry axis of the Kerr gravitational field) is the 
vanishing of the parameter L, i.e . L = 0. Assuming $ = 0, in equation (17), we 
can transform it into the Weierstrafi form of an elliptic curve by the following 
substitution 



-a/4 

Thus, we obtain the integral equation 

1 r aP' 



(19) 



(20) 



and this orbit integral can be inverted by the Weierstrafi modular Jacobi form 

5 

e = p(0M) (21) 

where A := — ^ — ,P' = {r^ + a^) and the Weierstrafi invariants take 
the form 

3.1 Exact solution for spherical polar null geodesies with 
a vanishing cosmological constant 

In terms of the original variables, the exact solution for the polar orbit of the 
photon ($ = 0) takes the form 

P(<^ + e) = ^cos2^-^(a"+/3") (23) 

where a" := a'/A'^ = = (3" := (3'/A'^ = Q" = Q'/A'^ = 

^5^. Also A' is given by the expression 

GMr 

^ ■■= — ~^or<.. (24) 

T-^ I 1 2GMr ^ ' 

— T ~T~ J- — 5 — 7~ 

Equation (23) represents the first exact solution of a spherical polar photonic 

orbit assuming a zero cosmological constant, in closed analytic form, in terms 
of the Weierstrafi Jacobi modular form of weight 2. 

^For more information on the properties of the Weierstrafi function, the reader is referred 
to the monographs [24, 25] , and the appendix of [10] . 
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The Weierstrafi invariants are given by 



{a" + I3"f Q"a" 



92 



9'i 



12 4 
1 (-a2 + Q)2 Q 



+ ■ 



12 a4A'4 4a2A'4 



216 48 48 

= 432i4^[-^«^-^'^^2 + ^'^'2^ + ^2^] 

The sign of the discriminant A'^ (A*^ = 9% — STgrf) determines the roots of 
the elUptic curve: A^ > 0, corresponds to three real roots while for A"^ < 
two roots are complex conjugates and the third is real. In the degenerate case 
A'^ = 0, (where at least two roots coincide) the elliptic curve becomes singular 
and the solution is not given by modular functions. The analytic expressions for 
the three roots of the cubic, which can be obtained by applying the algorithm 
of Tartaglia and Cardano [26], are given by 



(g^ + 2Q) A^ 
~ 12a2r2(GM/c2)2 

(g^ - Q) A2 
~ 12a2r2 (07^/02)2 
(2a2 + Q) A2 
~ 12a2r2(GM/c2)2 

(25) 

Since we are assuming spherical orbits, there are two conditions from the 
vanishing of the polynomial R{r) and its first derivative ^ . Implementing these 
two conditions, expressions for the parameter $ and Carter's constant Q are 
obtained 



$ = ,Q = 

a 



-^GM ^ ^2^ _ 3GM^2 ^ ^3 ^3(_4^2GM ^ ^(^3GM ^ ^-,2-, 



$ = Q = 

a(^-r) a2(^-r)2 

However, only the second solution is physical [22]. 



(26) 



^ These orbits are unstable since 4-^ > 0. However, they represent interesting new 

possible types of motion in the Kerr spacetime. They represent a non-trivial generalisation of 
the unstable circular closed orbit (photonsphere) in the Schwajrzschild black hole. 
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parameters 

OGalactic = 0.52 

$ = 0, r = 2.87313, Q = 25.8829 


half-period 
LO = 0.331117 


predicted dragging 
A(/,GTR = 1.32447 
- 75.88° per revolution - 273192 ^^'^^ 

^ revolution 


aCalactic = 0.9939 
$ = 0, r = 2.42451, Q = 22.3842 


w = 0.784737 


A0^^«- = 3.13895 
= 179.8° per revolution = 647455, 



Table 1: Predictions for frame dragging from galactic black hole for a photonic 
spherical polar orbit, for Kerr parameter ocaiactic = 0.52^^™,^ acaiactic = 
0.9939-^^^2^ respectively. The values of the radii arc in units of GMbh/c^) 
while those of Carter's constant Q in units of (GMbh/c^)^- The period ratios, 
T, are 2.33615i, 1.88276i respectively. 



The two half-periods u and u' are given by the following Abelian integrals 
(for A'^ > 0) [27]: 



dt 



1 V4i^ - 92t ■ 



93 



dt 



93 



(27) 



The values of the WeierstralB function at the half-periods are the three roots of 
the cubic. For positive discriminant A'^ one half-period is real while the second 

is imaginary ^. The period ratio is defined as t = 

An alternative expression for the real half-period oj of the WeierstraB function 



is: u 



%/ei-e3 

function 1 + + 
Thus « 



1, |Y3f|)> where F(a,/3,7,x) is the hypergeometric 



«(a+l);3(/3-H) 2 , 
1.2.7(7-M) ^ 



W- - 1 ■ 



ja^ + QHa^ + i-^+ry)^ 2 ^ 2 ' 2 ' ' + Q 



(28) 



After a complete oscillation in latitude, the angle of longitude, which deter- 
mines the amount of dragging for the spherical photon polar orbit in the general 
theory of relativity (GTR), increases by 



A(/.^™ = 4a; (29) 
We can also integrate the first and the third equation in (9). Then we get 



/ 



cdt 



■ 



1 



1 



dz 



2^ + 



1 - ^ ysv^v^o^yrri^ 



f 

Jo 



a^(l — z)dz 



2^/g^/iv/ (^^^l-(^) 



'^We organize the roots as: ei > 62 > £3. 

*Yet another equivalent representation for oj is u; : 



an c^a' 
,/0 I ^ 2GMi 



.f(i. i 1 
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Thus we obtain the foUowing exact expression for t 



ct = -4 



1 



+ 1 



2^4 + 1- 2GMr 



r(i)r(|) p/i 1 ^ 

r(l) V2'2' ' Q 



+ 



a- r(i)r(f) 



1 1 „ -a' 
' 2'2' '"q" 



2v/Q r(2) 

Similarly, if we integrate the differential equations for t and 4> we obtain 



(30) 



cdt 



or 

where we used Eq.(23) and the faet that, J p{(t>)d4) = —({(p), where ((z) denotes 
the WeierstraB zeta function. Also £ denotes a constant of integration. 

Assuming, that the centre of the Milky Way is a black hole and that the 
structure of spacetime near the region Sgr A* , is described by the Kerr geometry 
as is indicated by observations Eq.(l), we determined the precise frame dragging 
(Lense-Thirring effect) of a null orbit with a spherical polar geometry. We 
repeated the analysis for a value of the Kerr parameter as high as acaiactic = 
0.9939. Such high values for the angular momentum of the black hole, have 
been recently reported from x-ray flare analysis of the galactic centre [28] . The 
results are displayed in table 1. Let us also mention that in [29] it has been 
argued that an upper bound of a is given by a = 0.99616. 

3.2 Null spherical polar geodesies with the cosmological 
constant, Lense-Thirring effect and Appell hypergeo- 
metric functions 

We now derive an exact expression for the amount of dragging for a photonic 
spherical polar orbit in the presence of the cosmological constant, thus general- 
izing the results of the previous section. After a complete oscillation in latitude, 
the angle of longitude A<p, which determines the amount of dragging for the 
spherical polar orbit, is given by 



^P^-g^sin^g 



= a + 



A 

aP' _ 
A 



aP' 
A 



(31) 
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A<^GT^ = -4 



zdz 

7W) 



Oil 



1 r(i)r(i) 1 1 

Tl^XTT '^2'2'^'^'■ 



Q ' 3 ' 



1 ^(i)^(^^^dl,i,2,- 



^2' 2' 



(33) 



where f{z) = z{l - z){Q + z{Qa^^ + E^a^)){l + ^zf and ai = ^ - 
i£5!(g±£!), ^ _i££!(^£!)a!A^ function /?, /?', 7, x, y) is the 

first of the four Appell's hypergeometric functions of two variables x, y [20] ^, 

(34) 



™=o;^o (7,m + n)(l,m)(l,n) 
which admits the following integral representation 

/' u"-\l-uy-'^-\l-ux)-''{l-uy)-'''du = £^^lIllZ^Fi(a, /3, /?', 7, :r, y) 
^0 J- (7) 

(35) 

The double series converges when |a;| < 1 and \y\ < 1. The above Euler 
integral representation is valid for Re(a) > and Re(7 — a) > 0. Also 
r(p) = a;^'~-^e~^da; denotes the gamma function. 

For a zero cosmological constant (A = 0,(3i = 0) and we obtain the correct 
limit, 

,GTR 



(36) 



where ai = — in the limit of a vanishing cosmological constant. 

We can also obtain an exact expression for time. After a quarter of an 
oscillation in latitude the time elapses as 



ct 



L 



^ ill + 5iz)dz 



= 7fe'^^M2'2'l'^'- 



Q 



5i TT^ (3 1 ^ , Qo2| + SV -g^A' 



+ 7l2^M 2'2'^'^' 



Q 



(37) 



^The expression (A, k) = A(A + 1) • • • (A + k — 1), and the symbol (A, 0) represents 1. 
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In the limit A = 0, 71 = - 5 +° ' +\,Si ^ smdEq. (30) is recovered. 
The conditions from the vanishing of the polynomial R and its first derivative 
result in equations which generalize (26) and are provided in Appendix C. 



4 Non-spherical polar null geodesies 

In this case, the relevant differential equation for the calculation of deviation 
angle of light from the rotating black hole (or rotating central mass) is the 
following 

dct> _ 2aGMr 1 

dr- c^A ^ ^""^^ 
where the quartic polynomial R{r) is given by the expression 

R = r'^r\a''-Q)+'^-^{a' + Q)-a^Q (39) 

Expressing the roots of A as r+ , r_ which arc the locations of the event horizons 
of the black hole, and using partial fractions we derive the expression 

# _ ^? + ^- 



dr {r-r+)^fR {r-r-)^fR 



+ 



(r - r+)^{r - a){r - (3){r - 7)(r - S) 



(r — r_)-\/ (r — a)(r — /3)(r — — 5) 
where are given by the equations 

^^^^ (40, 

Also the radii of the event horizons are located at 

In order to calculate the angle of deflection we need to integrate the above 
equation from the distance of closest approach (e.g. from the maximum pos- 
itive root of the quartic) to infinity. We denote the roots of the quartic by 
a, /3, 7,5, a > /? > 7 > (5. Thus A(^ = 2 . We organize all roots in 
ascending order of magnitude as follows, 

a^> av> ai> ap (42) 
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We have the correspondence Q/^+i = oc, = P, Oiix-i = r+ = Q;/a-2, a^x-s = 7, = &■ 
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where = a^+i,ai, = otp = and a-i = ap^i,i — 1,2,3 and we have 

that a^-i > > 0:^-3. By applying the transformation 

u>z — 1 



or equivalently 

' - a^+2\ (r - afjt+i\ 



(44) 



where 

u; := ^^i^^ (45) 

we can bring our integrals into the familiar integral representation of Lauri- 
cella's Fd and Appell's hypergeometric function Fi of three and two variables 
respectively. Indeed, we derive 



+ 



H+ v^(i^(i-«^^)yrr7^ 

1/- APuj\a^+^ - a^+2) zdz 



+ 



f 

Jo 



H+ 




z){l — k\z)^\ — fJ,^Z 


-^^a;(a^+i 




dz 


H- 




-z){\-k'lz)^\-P?z 






zdz 


H- 




z){\- Ktz)^\- yP-z 



(46) 



where the moduli /t j_ , ^ are 



(47) 



Also 

= \/a;(a^+i - a^+2)(a^+i - Q!^_i)-v/";'+i ~ - "(u-s (48) 
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and a^_i = r±. By defining a new variable z' := uz we can express the angle 
in terms of Lauricclla's hypergeometric function Fd, 



A0°™ = 2 



r(3/2)r(i) 
r(5/2) 



A+V^(q:/.+i - Q;.+2) ^ /3 1 1 5 ^ 2 , 
^ if+ ^^1^2'2' '2'2'a;'''+'^ 

-2AgV^(a^+i - a^+2) ^ /I 1 1 3 1 2/ / 
+ ^= ^^U'2'^'2'2'^:;'''-''' 

AgV;U(a^+i - a^+2) „ 1 1 5 1 ,2 ,2\ r(3/2)r(l) ' 
^ H- ''y2'2' J r(5/2) 



(49) 

where the variables of the function Fd are given in terms of the roots of the 
quartic and the radii of the event horizons by the expressions 



1 


oin - 






6 


-/3 






a/x - 






5 


— a 








-< 


-1 




0- 


r± 






-< 


-1 




a — 


r± 


1^" = 




- (Xn- 

- ct^- 


-3 
-3 




P- 
a — 


7 
7 



(50) 



An equivalent expression is as follows 



GTR 



= 2 



-2AP^{a^+i - a^+2) A i 1 i 3 1 
H+ V2'2' '2'2'a; 



+ 



+ 



A^V^(a^+i-a^+2)/_^ /^i i i 3 1 

'•'^ ' 2' 2' 2' 2' w"" 



1113 1 

2' 2' '2' 2'^'"+ 



-2^gV^(a^+i - a^+2) ^ /I 1 1 3 1 2 ,,2\ 
+ ^= ^^V2'2'^'2'2';:'''-'^ ) 



+ 



H- 

1113 1 



1113 1 
2' 2' 2' 2' w' 



— T^r-^i I 7^> 1 2 



::^^-i2'2'^'2'^'^-''^"'2 



(51) 
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In going from Eq.(49) to Eq.(51) we used the identity which is proven in Ap- 
pendix C 

3 1 1 i ^ 1 .2 ,A r(3/2) 



2'2' '2'2'w' ; r(5/2) 

J_p /"l 1 1 ^ 1 .A r(V2) 

V2'2'2'2'a;'^ J T{3/2) 

1 ^ ^1 1 , 1 3 1 ,,2 A r(i/2) 



+ 32^i5 -,-,1, 



V2 2' ' 2' 2' ^"-^ y r(3/2) 



(52) 



The phenomenological apphcations of Eq.(49) for gravitational bending and 
lensing studies from a galactic black hole, as well as its generalization in the 
presence; of the cosniological constant will be the subject of detailed investigation 
in a future publication. 

4.1 Exact solution of timelike spherical poleir orbits with 
a cosmological constant 

Using the second and the fifth line of Eq. (4) , for L = and assuming a constant 
value for r, we obtain 

d(j) _ H a; 



d6 y/W 
-E^Ea + BAe 

where B :— °" '''a''"" ■ 

Similarly, using the third and fifth line we obtain 



(53) 



cdt S2(r2 + a2)2£ a^S^Ssin^ 



de ArVO' AeVO' 

TAg-a'^E'^Esm^e 



AeVW 



(54) 



andr:= ^^^-X ) 

Now using the variable z = cos^ 6, we obtain the following system of integral 
equations: 

r 772„ I? /o r Rn J- f 

^ ~' dz 



= f^dz+f. 
J y^} J 

J y^) ' J 



Ct = / ^ ' '-' dz - / ^ ; ' (55) 
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or 



/ 
/ 



(q'i + 0iz)dz 



ct = 



V7M 

{ji+diz)dz 



(56) 



where f{z) = z{l-z){Q + z{Qa^+E^a''E^-li^a^) + z\-n^a^^)){l + a^zf. 



Equation (55) is a system of equations of Abelian integrals, whose inversion 
in principle, involves genus-2 Abclian-Siegelsche modular functions. Indeed, this 
system is a particular case of Jacobi's inversion problem of hyperelliptic Abelian 
integrals of genus 2 [43]- [49] (see Appendix A for details). Then, one can express 
2; as a single valued genus two Abelian tliota function with argumcncnts t, 4>. 
However, since the polynonial f{z) of sixth degree posses a double root it may 
well be that the Abelian genus-2 theta function degenerates and the final result 
can be expressed in terms of genus- 1 modular functions. 

5 Frame dragging in spherical polar timelike geodesies 
with a cosmological constant 

After a quarter of oscillation in latitude the change of longitude is 



Also we have defined 



7i 



" ~2~ 2 A^^ 
l aE^{r^+a^)Ea^A 
2 3 

!^ i 1 

2 2 

a^AE^r^+a^E a^E^E 

6 Ar 2 



(57) 




Zi = - 



-Sa^E^E^ + 3a2 - a^QA - sJl2QKa^ + {^o?E-^E^ - Sa^ + o?QKf 



2a^A 



Z2 = 



-3a2£;2s3 + 3a2 - a^QA + ^ UQAa^ + {2,a?E'^E^ - 2.a? + a^QAf 



2a'^A 
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6 Spherical non-polar null geodesies 



In this section, wc investigate splierieal null geodesies with a nonzero value of 
the parameter $. Now one has to calculate the integral 



2 7 ^az^-z'^{a + p) + Qz 2] (^i - z)^az^ - z^a + p) + Qz 

(58) 

where 

a := -a^/3 := Q + $^ P = r'^ + - (59) 
The first integral can be brought into the Weierstrai3 form with the invariants 

92,93 

1 0"rr" 

_1_ [2{-a' + Q + $2)3 - 9Qa* + 9Qa'{Q + ^')] 

*_ _ GMr 

^ = r^'l, ^Mr (60) 

Now let us discuss the second integral in Eq.(58). We define a new variable 

u'lZ = V (61) 



53 

and 



then 



2 70 



dz $ 



2 70 {1 - z)^yaz^ - z^{a + p) + Qz 
$ dv 



Jo 



^ r(i)r(i) /I I 
2^ r(i) ^H^2' '2' 



2^2 

■^^^-F,(i,l,i,l,x.X2 



where ul- = i = 1, 2 and 



-a2 + Q + $2 _ w'4a2Q+(a2-Q-$2)2 

= 2a2 

-a^ + Q + $2 + ^4a2Q+(a2_Q_$2)2 

= 2^^ 
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parameters 








predicted dragging 




$ = 


= 1, Q = 22.693, r = 


= 2.7452 




-- 7.72736 


442.7° per revolution = 


1.59 X 10*^' T^ 

revolution 


$ = 


-1, Q = 26.9984, r 


= 2.99523 




-5.0551 = - 


-289.6° per revolution = 


1.04 X 10*^ ""'"l 

revolution 


$ = 


-3,Q = 23.0508, r 


= 3.2239 




-5.2014 = 


—298° per revolution = 


-1.07 X lO'^— 



Table 2: Predictions for frame dragging from a galactic black hole, with Kerr pa- 
rameter ttGalactic = 0.52 '^^'^" , for different values of photon angular momentum 
and Carter's Constant. The values of the radii and $ are in units of GMqh/c^, 
while Carter's constant Q in units of (GMbh/c^)^- 





parameters 








predicted dragging 


$ = 


1,Q = 16.1443, r = 


= 2.02083 


= 10.7355 


= 615° per revolution = 2.2 x 10^^^^^ 


$ = 


-1,Q = 25.8865, r 


= 2.73783 




-3.73503 


= 214° per revolution - 770405 


$ = 


-3, Q = 25.8628, r 


= 3.23713 


A<^GT=^ = 


-4.32779 = 


= 247.9° per revolution = 892671^^^^" 

revolution 



Table 3: Predictions for frame dragging from a galactic black hole, with Kerr 
parameter acaiactic = 0.9939 ^^°" , for different values of photon angular mo- 
mentum and Carter's Constant. The values of the radii and $ are in units of 
GMbh/c^, while Carter's constant Q in units of (GMbh/c^)^- 



Also 

^-=^' ^''^ 

Thus, we expressed the above integral in terms of Appell's first hypergeomctric 
function of two variables, ^1(0;, /?, /?', 7, a;, y). 

After a complete oscillation in latitude, the angle of longitude, which deter- 
mines the amount of dragging for the spherical non-polar photonic orbit in the 
General Theory of Relativity (GTR), is given by 

. PTR TT \^P-a\^/l 1 $ „ /I 1 



A V a c^a^ I F I _ _ 1 ' ^ I I A tt/? [ _ 1 _ 1 v v 

2x/Qv^^ + l-2^ U'2' 't.'J+ IsfQ^^ 'U' '2' 



(63) 



Orbits with A<^ > are called prograde and those with A^ < are 
called retrogade. 

As before we can obtain an exact expression for time. After a quarter of 
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oscillation in latitude 



dz 



2Vi ^az'^ - z{a + P) + Q 



+ 



dz 



a[a{l - z) - <1>] 



az^ 

1^2- 
2' 



+ 



A 

a' r(i)r(2) 
r(|) 



-z(a + /3) + Q 
— [(r^ +a^) -a$] -a- 



1 1 



VQ 



F 



i 1 1 !i 



P , 1 1 1 5 , , 
^i(2'2'2'2'"i'"^ 



(64) 



The above equation has the correct limit for $ = and reproduces the corre- 
sponding exact expressiion, Eqn.(30) for spherical null polar geodesies. Indeed, 

for $ = 0, ^ = ^g-, u'l = 1, and the Appell function has the following limit: 



r(i)r(2) 



^1 1 2' 2' 2' 2' ' Q 



r(^)r(i)^ 



r(2) 



1 1 -a2\ , , 
2' 2' 2'^) (65) 



Assuming that the centre of the Milky Way is a rotating black hole and that the 
structure of spacetime near the region Sgr A* is described by the Kerr geometry, 
we determined the precise frame dragging of a null orbit with a spherical non- 
polar geometry. The results are displayed in tables 2 and 3. 



7 Frame dragging in spherical non-polar null geodesies 
with a cosmological constant 

Including the contribution from the cosmological constant, the relevant differ- 
ential equation is 

where B' := S^aP. 

Introducing the variable z, we find that the angle of longitude that mea- 
sures the frame-dragging of a spherical non-polar photon orbit, after a complete 
oscillation in latitude is given by the exact expression 
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GTR 



= -4 



1$S2 1 vr /I , 1 , -«^A 1 u',' 



+ 



2 72^^1 V2' ' '2' ' 3u[ 'u'/u'i 
1 „ /I 1 a^A u'n 



ttFi -,1,-,1, 



2 

2^Q< ''^ V2'''2''' 3<'< 



B' 1 TT „ /I 1 W'n 



A, ^0^2 V2 2' 'u'l 



where = i = l, 2 and the roots ut are 

AE-«^ - 2AE-<]>r/-' + :^E-r/- + AE-<I>-r/- + QA(r - ;-iE-$- - 

Ml = 



(67) 



■W2 



and 



2a2 ((Aa2 - 2A$a + A$2 + 3) 52 + QA) 

AS2a4 - 2AS2$a3 + 3S2a2 + A52$2„2 ^ g^„2 _ 322^2 _ gg _ ^ 



2a2 ((Aa2 - 2A$a + A$2 + 3) 52 + QA) 



//a := 12Q((Ao2-2A$a + A$2 + 3)s2 + QA)a2 

+ ((a - $) (Aa^ - A$a2 + 3a + 3$) 5^ + Q (a^A - 3))^ 



(68) 



(69) 



Equation (67) involves three hypergeometric functions, Gaufi's F, Appell's 
Fi and LauriccUa's Fd- For vanishing cosmological constant the above expres- 
sion reduces exactly to (63). 

Similarly the time elapses after a complete oscillation in latitude is 



ct = -4 



lE2(r2+a2)[r2 + a2-a$] 1 



Ar 



-.F 



1 1 . w'o 



VQK V2'2' u[ 



1 0*52 fi 1 a2A u' 



--Fi 



.1 



+ 



2 VS< V2' '2' ' 3u[u[ 

a"S2 A _i 1 i 1 1 ^ !l 
' '2'''^'3ui'ni 



(70) 



8 Spherical timelike geodesies with L 

The relevant equation for integration is [11] 

1 f/- (aP L 

ad) = —- — X < — aE + 

2^z3a-z^{a + /3)+Qz {A 1-z 



(71) 
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where 



and P is provided from Eq.(lO) 
Let us define: 

n := 



a^a^{l~E^), (3 = Q + L^ 

dz j^ll J 



thusf = p(n+e), and ^" f , - , , 

— i — a-E) and e is a constant of integration. Now 



(1 - z)^/z3a-z^{a + P) + Qz 



(72) 
(73) 

(74) 



under the substitution (19) becomes 

La 



La 



La r 



di 



{w - i) V'4f - .92$ - .93 

p'(n)dn 



La 



{w - p(n + e)) V4p3(n) - 32 p(n) - 

{w - p(n + e)) 

CT(n + e-i;o) J 1 



Log 



a(n + e + t;o) 



+ 2n C(t;o) 



P'(^'o) 



(75) 



where w := ^{l — = p(i'o)- Also a{z) denotes the Weierstrai3 sigma 

function. Thus the equation for the orbit is given by 



= j d(j) = A"n - 



La 



Log ° + 2nC 

fT(n + e + vo) 



P'(^'o) 



(76) 



and p^'ivo) = 4p3(wo) - ff2p(wo) - ffa = 4'u;3 - - 93. In terms of the 
integration constants, w is given by the expression: 



w = 



(77) 



4 12 
Similarly using the first and third line of Eq.(9), we obtain for t the expression 



ct 



+ n an , 



n 



3 a 



+ — -C(n) (78) 



^The extreme black hole a = 1 solutions for L ^ have been investigated in [38] 
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An alternative exact expression for time t in terms of ordinary hypcrgcomct- 
ric of one- variable F and Appell's first hypergeometric function of two variables 
is as follows: 



ct 



-P + aL 



I o' o' ' 



2\/^ V2'2' '«i 



a^ii; r(i)r(i) /i i i u',\ 
+ m ^^U'" '2'''^4xj ^^^^ 

The variables are given in terms of the constants of integration and the 
Kerr parameter a by the expressions 



ui = 



a?{l -E^)+L^ + Q- V(-a2(l - E^) - - QY - Aa^jl - E^)Q 

2a2(l - £;2) 



U2 



a^jl -E^)+L^ + Q + vV^l - E^) - - - 4a2(l - £;2)Q 

2a?{l - E-^) 

(80) 

and u'i = l/ui,i = 1, 2. 

Similarly an alternative expression for the amount of dragging for timelike 
non-polar spherical orbits is as follows, 



A J,GTR _ /| V a ' c^a^ / " p ( ^ ^ 1 12. 



1 „ /-L I 2GMEr \ ^ /-, 



(81) 



where Vi = Yo = 

We have calculated the amount of frame dragging for galactic black holes for 
the values of Kerr parameter given in equation (1) and for fixed values for radii 
and Carter's constant Q. The invariant parameters L, E are determined by the 
two conditions for spherical orbits. We repeated the analysis for different values 
of the Kerr parameter [28] . The results are displayed in table 4. 



8.1 General solution for non-spherical geodesies in Kerr 
metric 

In the general case with non-zero cosmological constant, of non-spherical orbits 
one has to solve the differential equations 
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parameters 

a = 0.52, L = -2.03566, E = 0.957665, = 11 
a = 0.9939, L = -2.25773, E = 0.959284, Q = 11 
a = 0.99616, L = -2.25883, E = 0.959292, 0=11 



A(t>' 



GTR 
GTR 



predicted dragging 



GTR 



-6.06933 = 
-5.86166 = 
-5.86161 = - 



-1.2519 X 10*^' 
-1.2090 X 10*^ 

revolution 

■1.20904 X io6_^ics_ 



revolution 
arcs 



rBvnliitinTi 



Tabic 4: Predictions for frame dragging from a galactic black hole, with the in- 
dicated values of Kerr parameter acaiactic, for different values of test particle's 
angular momentum, a particular value for Carter's constant and for a fixed ra- 
dius r = 10. The values of the radii a,L are in units of GMbh/c^ , while Carter's 
constant Q in units of (GMbh/c^)^- The period ratios, r, are 3.0581,2.697431, 
2.696251 respectively. 



where R'{r) is a quartic polynomial for null geodesies (yu = in Eq.(5)) which 
is given by 

R' = E^ {S^ [(r^ + a^) - a$]' - A, [E\^ - af + Q] } (83) 
and 6' (6*) is given by 

& = E^\[Q + i^- am Ae - (54) 

sm J 

Note that this equation is an equation between two elliptic integrals, in the 
general case when all roots are distinct. The left hand side can be transformed 
into an elliptic integral with variable z or ^ in WeierstraB normal form, see 
Eq.(19). In order to solve this differential equation and determine r as a func- 
tion of 6, one can employ the theory of Abel for the transformation of elliptic 
functions [32], which was first applied in [11] for the case of timelike orbits. A 
detailed exposition of the theory of transformation of elliptic functions based 
on [32] can also be found in [11]. The interesting connection with modular 
equations [11] is outlined in Appendix B. 

We note at this point that the corresponding relationship for non-spherical 
timelike orbits in the presence of the cosmological constant relates a hyperelliptic 
integral (the polynomial R' is a scctic in this case) to an elliptic integral, and 
therefore involves the transformation theory of hyperelliptic functions. This fact 
was first observed in [11] and will be a subject of another publication. 



8.2 Transforming the geodesic elliptic integrals into Abel's 
form 



The transformation 



x^es + ^^l^ (85) 
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transforms the elliptic integral in Weierstrafi form into Abel's and Jacobi's form: 

f dx 1 f dx 

J ^y4{x-ei){x-e2){x-e^^~ Vei-~^J ^/(T^x^){l - k^x^) 

with A;2 = ^a^. 

61—63 

Similarly, the quartic can be brought to Jacobi's form y"^ = {1 — x^){l — kfx^) . 



Indeed, as we saw in section 3 we have 

d6 1 dz 1 d^ 



(87) 



where 

(88) 

In this case the three roots of the cubic 4^^ — 32^^53 ^'^e given by the expressions 
ei = (^-a^ + $2 ^ g + 3^0* + 2a2Q + Q2 _ 2a2$2 + 2Q$2 + $4^ 

63 = ^ (-a^ + $2 + Q _ 3v'a'' + 2a2Q + Q2 - 2a2$2 + 2Q$2 $4^ 

(89) 

Thus the Jacobi modulus fc^ = is given by 

ei— ea ° 



2Va^ + 2a2(Q - $2) + (q + $2)2 ^ ^ 

It has the correct limit for photon spherical geodesies with $ = 



Also 

1 2 



\/ei - 63 (a4 + 2a2(Q-$2) + (g + $2)2)i/4 
and it also has the correct limit for spherical photon orbits with zero $ 

1 2 



(92) 



(93) 
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Thus, we get 



r de 
J 7i 



(94) 



Applying Luchterhand's transformation formula [36] on the radial integral 
the Jacobi's form can be recovered: 

, , = (95) 

M^{l-x^){l-k'tx^) ^/{y-a){y-(3){y-^){y-5) 

where the Jacobi modulus fci and the coefficient M, are given in terms of the 
roots a, (3,^,5 of the quartic, by the following expressions 

^^ = ^^^H' M=^ia-,)iP-S) (96) 
^/{a-'y){|3-S) 

and the integration variables are related by 

1 - a: ^ - 5){y - a){y - (3) 
1 + x {a - p){y - ^){y - S) 



(97) 



Also wc assume that the roots a, f3, 7, S of the quartic are real and are organized 
in the following ascending order of magnitude: a > /3 > 7 > 5. 

We can also provide a nice formula for r in terms of the Jacobi's sinus 
amplitudinus function 



- 6){r - a){r - 0) _ ^ " (^/tI'^Q 
(a-/3)(r-7)(r-<5) 1 + gn (m / ^, /fci) 



(98) 



8.3 Equatorial geodesies including the contribution of the 
cosmological constant 

The equatorial geodesies (i.e. 6 = 'it/2,Q = 0), with a nonzero cosmological 
constant, may be obtained by Eq.(4) for the particular values of Q,0. The 
characteristic function in this case, has the form [11] 

W = -Ect+ f ^dr + L(l) (99) 
and the geodesies are given by the expressions: 



2wg wrieht f — = f rf^ = f 9x _ r _de_ 

^ ■> VR ■> ^(r-a){r-g)(r-f){r-S) ^ (l-:i2)(i_fc2^2) J s/B' 
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dr dX 



r 



.2 



_ a{l + \a'mE{r^+a')-La) 1 
"rfA - (i-Ar-2)(,2 + «2)_2^ +(^ ai;)(l + -aA) 

cr- = + {l + -aK)a{L-aE) 

(100) 



where 



i?' = (1 + ia^A)^ [((r^ + a'')E - ahf - A,((L - aEf)] (101) 
o 

for null geodesies and 

R' = {l + \a'^Kf [((r^ + a'^)E - ahf - A^((i - aSf)] - Arifi^r'^) (102) 

for timelike geodesies. 

Equatorial orbits are of particular interest for various astrophysics appli- 
cations. The exact solution of circular equatorial orbits with a cosmological 
constant was presented in [11]. In what follows, we shall concentrate on the 
cases of non-circular equatorial timelike and null orbits that describe motion 
of test particles and photons under the assumption of a vanishing cosmological 
constant. We shall derive new exact expressions for the precession of perihelion 
or periapsis for the orbit of test particle in the gravitational field of Kerr, as 
well as for the deflection angle of a light ray from the Kerr gravitational field, 
in terms of multivariable hypergeometric functions. In the latter case, we apply 
the exact calculation obtained for determining the bending angle of a light ray 
from the gravitational field of the galactic centre of Milky- Way assuming that 
the Sgr A* region is a supermassive rotating black hole for various values of the 
Kerr parameter which are supported by recent observations and of the impact 
factor. The more general case, in the presence of the cosmological constant is a 
task for a future publication. 



8.4 Exact solution of timelike equatorial geodesies 

We now proceed to determine the exact expression for the precession of equa- 
torial timelike orbits in Kerr spacetime. We have r'^{r) = y/R. This can be 
rewritten as 

= + + (103) 



By defining a new variable u = l/r we get the following expression 

^5— (L - aEYu'^ - (1 + a^u : 



u-\^ = E^ + a^E^"" -L\'' + =^ {L - aEfu^ - (1 + a^u^ _ ^^u) = B* {u) 

(104) 
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Similarly = w^^^M where 



1GM. 

A{u) = L + uas{aE - L), D{u) = I + a'^u'^ - asu, as := — ^ (105) 



Thus, we obtain the differential equation 

dcj) A{u) 1 



du D{u) ./Wiu) 

We now write 



(106) 



D(u) u+ — u U- — u 
where «+ = , = -t- and 



r± = — ±J(^— j -a^ (108) 



where the quantities A+, A_ are given by 



U- — u+ 
U- — u+ 



(109) 



For the calculation of the perihelion (periapsis) precession of a test particle 
in orbit around a rotating mass we need to calculate the integral A^*^™ = 

2j^^d<j). Then 

f dd> = I '^'"'^^ . ^ . ^ 



du'A_ 1 1 



(110) 

and we have defined u' := u^-. Using new variables 
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wc can bring our expression to the integral representation of Appell's first hy- 
pergeometric function 

. .GTR ^ _^ 1 \ A+ /^i T i 1 "2 - »3 r(l/2)^ 



GMr- , " ^ 1 o' o' GMr_ ) ~> 



Wc can now provide an exact expression for time t. The first differential 
equation in Eq.(9) for equatorial orbits can be written in terms of the variable 
u as follows 

E{l + a''u'')+aasu\aE-L) 

Subsequently by dividing with ii we derive the equation 

cdt _ E{1 + a-tr) + (insir''{(i.E - L) 

E Ea^ auas{aE — L) 

v?D{u)s/B^ D{u)s/Wiu) D{u)s/BH^ 

(113) 

By integrating we can express t in terms of Appell's and Lauricella generalized 
hypergeometric functions 



(111) 



V (^) 

+ -nnZ r-f*! 7^)1)7^>1) 



A'_^_ 1 1 J 1 p ( 1 r, 1 1 1 ^2 ~ '"3 '"3 ~ '"2 ^ ~ f*2 



,' n,' nJ 



-F^ 2, 1, 1, ^Ji_ \ , ^ (114) 



GMr- „/ -^19' ' '9' ' „' 'GMr- 

8.5 Exact solution of null equatorial geodesies 

In this case, we arrive at the differential equation first derived in [14] 

_ $ + uas{a - $) 1 



(115) 
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where the cubic polynomial (u) is given by the expression 

B^{u) = as{^ - afu^ + u^{a^ - + 1 (116) 
In order to calculate the angle of deflection is necessary to calculate the integral: 

Using partial fractions as in the timelike case we obtain an elegant exact 
expression for A^'^"'"'^ in terms of Lauricella's fourth, hypergeometric function 
of three variables Fjj, and Appell's first hypergeometric function of two variables 



A(l)' 



GTR 



+ 
+ 



^-'s 'V' '2' '^-u','<-u',) r(i) _ 



1 



u[ — Mg I as if -a)^ 



~^+\ ~ S"T7T77 x-2-Fd ( TTi 1, 7:, 7:, 7:5 ■ 



V" 



-3 ^^-4) V 
A I 1 op 1 1 1 3 

~ ^-1 77777: X^^^ o'l'o'o'o' 



(117) 



or 



1 _^ 1 1 3 -u'-i W 




2' '2'2'2'g^_v'<-«^'«'2-«'3 



3 

I p I 1 1 1 1 "9 - m'. 



o' 9' GMr- I \i' — q,' 



3 . 



— Uo „ / 1 , 1 1 3 —Mo —U'r, 



/ "3 p \ L ^ L t ■! ~ "3 

\/ ^ o' ' o' o' o' GMr_ / ' „./ 



2' '2'2'2'^_^'«;-«-«',-«^ 



(118) 
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parameters 




A^GT^ - 


predicted deflection 


aCalactic = 0.52, $ = 


5 


TT = 1.84869 = 105.9° = 


= 381319arcs 


aCalactic = 0.52,$ = 


10 


A</)«™ - 


- TT = 0.537 = 30.77° = 


110792arcs 


aCalactic = 0.52, $ 


15 


A(/)«™ - 


- TT = 0.320 = 18.34° = 


66015. Bares 


■ = 0.52, $ = 


20 


A0«^^ 


-TT = 0.228 = 13.08° = 


= 47093arcs 


aCalactic = 0.52,$ = 


40 


A0«^^ 


- TT = 0.1065 = 6.10° = 


= 21973arcs 



Table 5: Predictions for light deflectiom from a galactic rotating black hole with 
Kerr parameter acaiactic = 0.52 '^^'^" . The values of the impact parameter $ 
are in units of ^^^^ . 



During the derivation of Eq.(118) we have used at an intermediate step of the 
calculation the identity 

1 / 1 1 3 m'^ u'^\ 

z-tn \ i, i. 



1 „ / 1 1 1 3 



(119) 

The quantities A± are now given in terms of the impact parameter $ and 
the Kerr parameter a by the expressions 

A± = ^ '"^ (120) 

The angle of deflection 6 of light-rays from the gravitational field of a Galactic 
rotating black hole or a massive star is defined to be the deviation of Acj)'^'^^ 
from the transcendental number tt 

6 = A<^GTR _ ^ (121) 

We have calculated the deflection angle S of light rays from the gravitational 
field of a galactic rotating black hole for different values of the Kerr parameter a 
and the impact parameter $. The results are displayed in figure (1) and tables 
5,6. It is clear from figure (1) that especially for smaller values of the impact 
parameter $, there is a strong dependence of the deflection angle on the spin of 
the black hole. This has implications for gravitational lensing studies and can 
lead in principle to an independent measurement of the Kerr parameter at the 
strong field regime. 

The roots of the cubic are organized as u'^ > ^2 > > U3. 
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Figure 1: Deflection angle S for various values of impact parameter $ and for 
two values of the Kerr parameter a for the rotating black hole. 



32 



parameters 

aoaiactic = 0.9939, $ = 5 
acaiactic = 0.9939, $ = 10 
acaiactic = 0.9939, $ = 15 
aCaiactic = 0.9939, $ = 20 
acaiactic = 0.9939, $ = 40 



predicted deflection 
A0°™ - TT = 1.317 = 75.47° = 271710arcs 
A0°™ - TT = 0.497 = 28.49° = 102581arcs 
A(j)'^™ - TT = 0.306 = 17.57° = 63254arcs 
A0°™ - TT = 0.2217 = 12.70° = 45728arcs 

GTR 



Act) 



TT = 0.1052 = 6.026° = 21694.3arcs 



Table 6: Predictions for light deflectiom from a galactic rotating black hole with 
Kerr parameter acaiactic = 0.9939^^^. The values of the impact parameter $ 
are in units of '^^^^ . 



8.5.1 Null equatorial geodesies with a cosmological constant 

In this case, the generalization of equation (115) is given by 



# _ asv?{^ -«) + ($- a)(-4)(l + o?v?) 



du [(^2 - f ) (1 + a2u2) _ asu^] ^/W^ 

where B^{u) is still a cubic polynomial 



(122) 



B^iu) = 1 + a^u" - + _ af + ^(1 + aV)($ - af (123) 

For a vanishing cosmological constant the above cubic polynomial reduces to 
Eq.(116). 



9 Conclusions 

In this work, we have investigated the motion of a test particle and light in the 
gravitational field of Kerr spacetime with and without the cosmological con- 
stant. Wc have derived a number of useful analytical expressions for measurable 
physical quantities. 

In the case of null orbits, we solved exactly the geodesic equations for spher- 
ical polar and non-polar photon orbits. The exact solution for the orbit of a 
photon with zero angular momentum $ and vanishing cosmological constant 
was provided by the Weierstrafi elliptic function p{z). 

The exact expressions that determine the amount of frame-dragging (Lcnsc- 
Thirring effect) for the corresponding photon orbits, assuming vanishing cosmo- 
logical constant, were written in terms of Weierstrafi function real half-period or 
equivalcntly in terms of Gaufi hypergeometric function F for a photon spherical 
orbit with $ = and Gaufi hypergeometric function F and Appell's generalized 
hypergeometric function of two variables Fx for photon orbits with $ ^ and 
constant radius. The corresponding expressions in the presence of the cosmolog- 
ical constant were given in terms of Appell's hypergeometric function Fy for the 
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case of spherical null polar orbits, and in terms of Lauricella's liypergeometric 
function of three variables Fd, Appell's Fi and Gaufi ordinary hypergeometric 
function F, in the case of spherical photonic orbits with nonvanishing value for 
the invariant parameter $. 

We subsequently applied our exact solutions for the determination of the 
Lense-Thirring effect that a photon experiences in a spherical polar and non- 
polar orbit around and close to our galactic centre, assuming the latter is a 
rotating black hole whose surrounding spacetime structure is described by the 
Kerr geometry as supported by recent observations. We repeated the analysis 
for various values of the Kerr parameter. 

We also solved exactly non-spherical polar null unbound orbits. We derived 
analytical results for the deflection angle of a light ray from the gravitational 
field of a rotating black hole's pole. The resulting expression for the deflection 
angle was written elegantly in terms of Lauricella's hypergeometric function Fd ■ 

We then investigated, non-circular orbits confined to the equatorial plane 
(timclikc and null-like equatorial geodesies) for which the vahic of Carter's con- 
stant invariant parameter vanishes. For the case of a vanishing cosmological 
constant, we derived an exact expression for the amount of relativistic preces- 
sion for a test particle in a timelike orbit, around a rotating central mass. The 
corresponding novel expression was given in terms of Appell's first hypergeo- 
metric function of two variables Fi. The application of this exact solution as 
well as of those that describe non-spherical orbits not necessarily confined to 
the equatorial plane, for the determination of the effect of rotation of central 
mass on the perihelion precession of a test particle (Mercury around Sun) or 
periapsis precession for a star such as S2 in a high eccentricity orbit around the 
galactic centre [42] is beyond the scope of this work and will be the subject of 
a future publication [51]. 

We have also derived an exact expression for the defiection angle of a light 
ray from the gravitational field of a rotating mass (the Kerr field). The cor- 
responding expression was given in terms of Lauricella's Fq and Appell's Fi 
generalized hypergeometric functions. We applied this calculation for the bend- 
ing of light from the gravitational field of the galactic centre of Milky Way, 
assuming the latter is a supermassive Kerr black hole for various values of the 
Kerr parameter and the impact factor. We emphasized in the main text the 
strong dependence of the bending angle on the Kerr parameter for small values 
of the impact factor. These results should be useful for gravitational lensing 
studies, where one treats the black hole as a gravitational "lens" [50], especially 
in the strong field region, when the bending angle can be very large. 

The synergy between theory and experiment for probing and measuring rel- 
ativistic effects is going to be one of the most exciting and fruitful scientific 
endeavours. 
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A Definitions of genus-2 theta functions that 
solve Jacobi's inversion problem 

Riemann's theta function [48] for genus g is defined as follows: 

m,--- ,ng 

where fln^ := Ounf + • • •2r2i2nin2 + ■ ■ ■ ;>ncl un := uiui + • • -UgUg. The 
symmetric g ^ g complex matrix whose imaginary part is positive definite is 
a member of the set called Siegel upper-half-space denoted as ■ It is clearly 
the generalization of the ratio of half-periods r in the genus 5=1 case. For 
genus g = 2 the Riemann theta function can be written in matrix form: 



E 



7ri( m ^2 ) 



( f^ll 


^12 ] 


(:) 


I ^^12 


^22 ] 





+2fri( m 712 



Ui 
U2 



e 

n\,n2 

(125) 

Riemann's theta function with characteristics is defined by: 

Q{u\q,q') := g2niu{n+q')+iT'^{n+g'f+'^-^iq(.ri+q') (126) 

herein q denotes the set of g quantities qi, • • • ,qg and q' denotes the set of g 
quantities q'l, - ■ ■ , q'g- Eq.(126) can be rewritten in a suggestive matrix form: 







q' 
q 



(w,0) = ^ e'''*("+«')"("+«'^+2'^'*("+«'^("+«\ q,q'€Q^ (127) 



neZ9 



The theta functions whose quotients provide a solution to Abel-Jacobi's 
inversion problem are defined as follows [49]: 

6{u; q, q') := gau'^+2hu(n+q')+b{n+g')^+2iTrg{n+g') (128) 

where the summation extends to all positive and negative integer values of 

the g integers ni, • • • ,ng, a is any symmetrical matrix whatever of g rows and 
columns, h is any matrix whatever of g rows and columns, in general not sym- 
metrical, b is any symmetrical matrix whatever of g rows and columns, such 
that the real part of the quadratic form bm? is necessarily negative for all real 
values of the quantities mi,-- - ,mg, other than zero, and q,q' constitute the 
characteristics of the function. The matrix b depends on ^g{g + l) independent 
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constants; if wc put iirfl = b and denote the ^-quantities hu by inU, we obtain 
the relation with Riemann's thcta function: 

0{u;q,q')=e''^"e{U;q,q') (129) 

The dependence of genus-2 theta functions on two complex variables is de- 
noted by: 0{u;q,q') = 9{ui,U2;q,q'), the dependence on the Siegel moduli 
matrix O by: 6{ui,U2,0,;q,q'). To every half-period one can associate a set of 

characteristics. For instance, the period u"'"! ~ 5 ^ | g 

theta function of two variables with zero characteristics, i.e. 9{u) = 9{u; 0, 0) = 

P Q (u, Ct). Also, Weierstrafi had associated a symbol for each of the six 

odd theta functions with characteristics and the ten even theta functions of 
genus two. For example, 9{u) is associated with the Weierstrafi symbol 5 or 
occasionaly the number appears as a subscript, i.e. 9{u)r,. 

Let the genus g Riemann hyperelliptic surface be described by the equation: 

= 4:{x — oi) • • • (a; — ag){x — c){x — ci) ■ ■ ■ {x — Cg) (130) 

For g = 2 the above hyperelliptic Riemann algebraic equation reduces to: 

y2 = 4(2- _ ai){x - a2){x - c){x - ci){x - C2) (131) 

where ai, 02, cci, C2 denote the finite branch points of the surface. 

The Jacobi's inversion problem involves finding the solutions, for Xi in terms 
of Ui, for the following system of equations of Abelian integrals [49]: 

Xl.ai , I ^g^O^g 

Wl + • • • + Ml = Ul 

: + ••• + : : 

where = T f , = /; ^, ' • • , = /; 

For g = 2 the above system of equations takes the form: 

dx dx 

/ — + / — = Wl 

J ai y J a2 y 

X dx f^^ X dx , , 

+ / = U2 133 

y Ja2 y 

where ui,U2 are arbitrary. The solution is given by the five equations [49] 

'2 , = A{b-xi){b-X2)---{b-Xg) 
9 (tt) 

= A{b-X,){b-X2) 



(b- Xi_)(b - X2) 
^e-^PP'f'{b) 
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where /(x) = {x—ai){x — a2){x — c){x — ci)(x — C2), and e'^*"^^ = ±1 accordingly 
as u'''" is an odd or even half-period. Also b denotes a finite branch point and 
the branch place a being at infinity [49]. The symbol 6{u\u^'°') denotes a genus 2 
thcta function with characteristics: 9{u:q,c/) [49], where ?i. = {ui,U2), denotes 
two independent variables. From any 2 of these equations, eq.(134), the upper 
integration bounds xi,X2 of the system of differential equations eq.(133) can 
be expressed as single valued functions of the arbitrary arguments Ui,U2. For 
instance, 

1 



xi = ai + 



Ai{x2-ai) e'^{u) 



(135) 



and 



X2 



[0.2 - ai)(a2 + ai) + 



6\u) 



2(ai - 02) 



± 



where 



(aa - ai)(a2 + ax) + 37— 



1 g'^(M|M°2.°) 



4(ai - 02)77 



2(ai - 02) 



02 6|2(m|u"1'") ai 6'2(w|u"2.«) 



77:^0201(01-02) , ^ , 2, 



^1 6»2(u) ^2 6»2(u) 



(136) 
(137) 



Also, Ai = ±- 



The solution can be reexpressed in terms of generalized WeierstralB functions: 



(^ 2) P2,2(w) ± y'pi,2W +4p2,l(w) 

fc' = n ' ^ = 1) 



(138) 



where 



P2,2(m) = 



(ai - a2)(a2 + ai) - ^ + ^ 



and 



P2,l(w) 



Ol — 02 



-0102(01 - 02; - -J- — -^rr:r\ r — — 



(«) 



(139) 



(140) 



oi — 02 

Thus, a;i,a;2, that solve Jacobi's inversion problem Eq.(133), are solutions of a 
quadratic equation [44, 49] 



Vx^ - U'x + U" = 



(141) 



where U, U' , U" arc functions of mi, 112- In the particular case that the coefficient 
of x^ in the quintic polynomial is equal to 4, [/ = 1, U' = p2,2(u), U" = p2,i('o). 
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The matrix elements hij , flij can be explicitly written in terms of the half- 
periods J7^^'". For clarity, U^*'''^ = Q xdx/y,Ul*'^^ = Qdx/y etc. The 
roots have been arranged in ascending order of magnitude and are denoted by 
623, e2g-i, • • • ,60,5 = 2, so that 621,621-1 are respectively, Cg_i+i, Og-j+i and 
60 is c. For instance, the matrix element /in = , (!4,e3..e2.4 — ■-.e2,eirre4,e3^ x ttj, 

Wmie Sill — jy*!2 .«1 [;-«4 .«3 _[/*'2 .«1 [j«4 >«3 • 

A.l A particular inversion problem of Jacobi 

An indefinite elliptic integral of the third kind can be regarded as a special form 
of a hyperelliptic integral 



L 



{a + (3x)dx ^^^2) 



^Jx{l - x){l - k'^x){1 - A\)(1 - n'^x) 

in the special case when two of the moduli are equal, e.g. X = fi, and therefore 
one can consider the following Jacobi's inversion problem 



{a + l3x)dx ^ {a + Px)dx 



r 

Jo (1 - A^a;)i/^(F^^)(1^^«^ -^0 {1 - X^x)^/x{^^^x)iT^^l^ 
ja' + (3'x)dx ja' + P'x)dx 

Jo (1 - X^x)y/x{l - X){1 - K?X) Jo {1- X'^X)^/X(T^^X}(T~^ 

(143) 

For a convenient choice of the constants a,l3,a',l3' the solution of the above 
Jacobi's inversion problem can be expressed in terms of genus 1, Jacobi theta 
functions 



±K^/XlX2 = 

4v(i-^i)(i-a=2) = 

k' 

^^/{l-K^Xl){l-K^X2) = 



e{a)e-'"e{u - a) - e''e{u + a) 



ei{a)e-"t 


h{u- 


a) + 




h{u + a) 




h{u- 


a)- 




hiu + a) 


6ii(a)e-"< 


h{u- 


a) + 




h{u + a) 




h{u - 


a)- 




?3{u + a) 




h{u- 


a) + 




h{u + a) 



(144) 



B Transformation theory of elliptic functions and 
modular equations 

One of the applications supplied by the transformation theory of Elliptic func- 
tions, which is of great importance in Number theory, are the modular equations 
described below [32, 33, 35]. 
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For a rational solution of the differential equation 



dy d. 



V(l-y2)(l-e2y2) ^(l_^2)(i_e2^2) 

the necessary conditions among the periods 

K{ei) = aoCK{e) + aiCiK'{e) 
iK'{ei) = boCK{e) + biCiK'{e) 



(146) 



with the period ratios (moduli) of the associated modular theta functions being 
given by 

r = ^^^±^ (147) 

are also sufficient, when 

Oo6i — aifeo = n (148) 
is a positive integer number. The integer n is called the degree of transformation. 

Equation (148) for ao,6i,ai,foo € Z when viewed as the determinant of a 
matrix e GL{2, Z), sometimes is called a modular correspondence of level n. 

It can be shown that the inequivalent reduced forms of modular correspon- 
dences 

ao ai 
bo bi 

are of the form 

q 

lee q' 

where q a positive part of n represents, := ^, and < £. < q' — I- For instance 
for n = p a prime number, there are p + 1 inequivalent reduced forms of the 

form ^"^ 

10\/10\/1 0\ / 1 \ f p 
p ) '\ 16 P ) '\ 16.2 p ) 16(39 - 1) P J \ o 1 

Also the multiplication factor C in Eq.(145) is given by 

which for a degree of transformation that is a prime number (n = p) is equal 
to or (1/n) ^^p-j ■ Equation (149) can be re-expressed in terms of Jacobi 

theta functions as follows: 



^^In a more familiar notation these classes of inequivalent reduced forms are 

a b 
d 



'^i ~ { n ^ ] ,ad = n, {a,b,d) = 1, < b < d and a,b,d E Z 
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The modular equations are equations relating the Jacobi modulus e(pr) to 
e(r) which are of the form 




(151) 



where denotes the Legendre symbol . Equivalently, modular equations 

can be written in terms of the absolute modular invariant function j(r), and 
relate the reduced absolute modular invariant j* to j by polynomial equations 
of the form 

^p(r,i) = (152) 

where j* ~ j-Ui = j (^^J^) The explicit form of ^2{j*,j) = 0, has been given 
in [41]. 



C Conditions for radii for spherical null geodesies 
with a cosmological constant and differential 
equations of Appell's function 

The conditions from the vanishing of the polynomial R and its first derivativere- 
sult in the equations which generalize (26) 

r^A^a^ + Ya^ +Xa± V3^/h 
~ r3A2a4 + (2A2r5 + 6^Ar^ - 9r + 9^)a2 + r^A{Ar^ - 3r + 6^) 

"(153) 

where 

/i := r2(Zo6 + r(3^/<: + r(2AV-15AV + 9Ar2 + 54)) 
+ r2(54(Ar2 + 2)(^)2 + UrK,^ + r'K,)a' 
+ 3(3^-r)r«A(Ar3-3r + 6^)) 
Indeed, for a vanishing cosmological constant it has the correct limit $ — > 

2, 2 

° derived in the previous subsection. 



Then the parameter Q is given by the expression 



Q 



(154) 



14 _2_' 
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where 



AV + 6^Ar*-9 



GM 



Y 

Z 
K 



2„5 



2A"r 



I 3^6 



GM 



(6Ar2 + 9) 



5AV - QAr^ - 36 

A^r^ - Vlt^r'^ + 18Ar^ + 27 

2A^r^ - 6Ar^ - 9 



(155) 



C.0.1 Proof of Identity Eq. (52) 

z'dz' 1 1 



dz' 



dz' 



1 - K'f^' 

1 



1 - k'^z' 



_ _ 1 /■^ d^' 1 1 /■^ rf^' 1 1 



(156) 



Picard, had developed a theory for finding solutions of the system of differen- 
tial equations that the Appell hypergeometric function obeys. More precisely he 
showed, by direct substitution, that solutions are provided by definite integrals 
of the form [30] 



J a 



(157) 



where g, h denotes two of the quantities 0, 1, x, y, oo and we have the correspon- 
dence 

6i = 1 + /3 + /3' -7,62 = 7-",/" = 1-/3,A = 1-/?' (158) 

This is the generalisation of Kummer's work who found that the standard 
hypergeometric equation has 24 solutions [31]. The system of linear differential 
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equations of Appell's function Fi is 

d'^F dF 
x{l -x){x- y)-^ + [-fix -y)-{a + f3 + l)x^ + {a + f3-0' + l)xy + 0'y] 

dF 

-Py{l-y)-^-a(3{x-y)F, = 0, 



d F dF 
]-^ + [^{y-x)-ia + p' + l)y^ + {a + p' -(}+l)xy + ^x]-^ 



BF 

-f3'x(l - x) — ^ - af]'(y - x)Fi = 0, 
ox 

. d^F^_^,dF, ^dF, 

dxdy dx dy 



(159) 



For instance the following integral is represented as follows 



OO 

y/3+/3'-7(y _ l)7-«-l(y _ a.)-/3(y _ y)-P'du = 



B(l + /? + /?' - 7, 7 - a)x-^y-'^'Fi (l + p + p' - + p + p' - a,-,- 

V X y^ 

(160) 

and B{p, q) = ^^'^^^^^^ denotes the beta function. 
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